Practical scheme for a light-induced gauge field in an atomic Bose gas 
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We propose a scheme to generate an AbeUan gauge field in an atomic gas using two crossed laser 
beams. If the internal atomic state follows adiabatically the eigenstates of the atom- laser interaction, 
Berry's phase gives rise to a vector potential that can nucleate vortices in a Bose gas. The present 
scheme operates even for a large detuning with respect to the atomic resonance, making it applicable 
to alkali-metal atoms without significant heating due to spontaneous emission. We test the validity of 
the adiabatic approximation by integrating the set of coupled Gross-Pitaevskii equations associated 
with the various internal atomic states, and we show that the steady state of the interacting gas 
indeed exhibits a vortex lattice, as expected from the adiabatic gauge field. 
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A major motivation of current research with atomic 
quantum gases is the prospect of simulating the physics 
of condensed matter systems (see, for example, [1 ). One 
objective is to gain new insight into strongly correlated 
states of matter such as the fractional quantum Hall ef- 
fect exhibited by electrons in a magnetic field [2 . To 
study such systems with cold atoms, one needs to mimic 
a gauge field with which the neutral atoms interact as 
if they were charged. A well-known method consists in 
setting the gas into rotation and transforming to the co- 
rotating reference frame, where the physics is governed 
by the same Hamiltonian that describes charged parti- 
cles in a magnetic field. At the mean-field level the cor- 
responding gauge field leads to the nucleation of vortices, 
which has been observed experimentally (see [3] and ref- 
erences therein). However, technical difficulties have so 
far prevented the strongly correlated regime from being 
reached, and it is worthwhile to contemplate alternative 
approaches. 

A promising possibility is the use of geometric poten- 
tials induced by laser fields. To comprehend the under- 
lying mechanism, consider atoms with a ground state ^, 
which is degenerate in the absence of any external per- 
turbation (angular momentum > 0). Neglecting in a 
first step the population of excited levels, the atomic en- 
ergy eigenstates |?/^q;) (a = 0, ...,2J^) in the presence 
of the laser fields are linear combinations of the Zee- 
man substates |^, J^, m) (m = — J^, . . . , Jg). Both the 
eigenenergies Ecy and the states |?/^ck) depend on position 
via the spatial variation of the light fields. If the centre- 
of-mass (COM) motion is slow enough, an atom initially 
prepared in one of these internal states, say IV^o), will 
remain in this state and its COM wave function will ac- 
quire a geometric phase [4]. This so-called Berry phase 
[5 corresponds to a gauge field {A, U} appearing in the 
effective Hamiltonian of the COM motion: 

g= ^"!y^' +^o(r-) + ^(r-), (1) 
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where M is the atomic mass, 

(2) 

and with the charge g = 1. Several configurations to 
create gauge fields with laser beams have been proposed 
and analysed [6 , 7 , 8 , 9] [lOl [IB [H [13l [M] , and experi- 
mental evidence of these potentials has been provided in 
the context of cold atomic gases [151 US]- particular, 
the proposal of [13 uses two counterpropagating beams 
driving the two transitions of an atom with a A level 
scheme [Fig. [TJa)]. Such a configuration emerges when 
the two beams are circularly polarised with positive and 
negative helicity, respectively, and both the ground {g) 
and the excited (e) states have unit angular momentum 
( = Je = 1) [17 . In this case |?/^o) can be chosen as 
the "dark" superposition state of |^, J^, m = ±1). Unfor- 
tunately, this scheme is not applicable to the commonly 
used alkali-metal atoms even though they may have some 
hyperfine states with the proper angular momenta. The 
hyperfine splitting of the relevant excited manifold being 
typically smaller than 100 linewidths, non-resonant cou- 
pling to other hyperfine states with Je 7^ 1 will destabilise 
the superposition state, and heating due to spontaneous 
emission will eventually mask the effect searched for. 

In this paper we first propose a simple planar scheme 
which overcomes the above mentioned handicap. It in- 
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FIG. 1: (a) A-type system used for a simple modelling of the 
atom-laser coupling, (b) Atomic scheme relevant for a 3/2 
nuclear spin, as for ^^Rb or ^^Na. 
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corporates two crossed laser beams whose frequencies are 
far-detuned from the atomic resonance transition, which 
reduces the rate of spontaneous emission processes to a 
neghgible value. Our scheme produces an effective gauge 
field that varies smoothly over the area of the atomic 
cloud and allows for the observation of several vortices 
in a cold interacting spinor Bose gas. In a second step, 
we determine the ground state of the gas by solving nu- 
merically the coupled set of Gross-Pitaevskii equations 
(GPEs) associated with the various internal (ground and 
excited) atomic states. We show that the ground state of 
the system indeed contains several vortices with the pre- 
dicted surface density. In contrast to previous approaches 
[18] where the gauge fields ([2| were explicitly included in 
the Hamiltonian, we do not a priori assume the adiabatic 
limit to be valid. Hence, our study can be considered as 
a validity check of this approximation for an interacting 
gas. In the last part of this paper we investigate whether 
the vortices can be nucleated by slowly turning on the 
geometric light potential. 

The necessity of using off-resonant lasers strongly lim- 
its the choice of schemes that lead to a nontrivial gauge 
field. To be explicit, we first remark that a frequency 
detuning large compared to the hyperfine splitting of the 
excited manifold causes the nuclear spin to become irrel- 
evant for the description of the excited state. Consider 
for example the Di transition between the states S1/2 
and Pi/2 of an alkali- metal atom, which is sketched in 
Fig. Ub) for the case of a nuclear spin 3/2, as for ^^Rb 
or ^^a. Since only the electronic angular momentum is 
relevant for evaluating the matrix elements of the atom- 
laser coupling, the energy level structure can be regarded 
as a simple Jg = 1/2 ^ Je = 1/2 system. At the same 
time this restricts the polarisations of the laser beams 
that can be applied to create the internal superposition 
state: A spin fiip from m = — 1/2 to m = +1/2, for 
instance, corresponds to an angular momentum trans- 
fer of so a two-photon transition requires light with 
circular (cr_ or cr+) as well as with linear (tt) polari- 
sation. Such a configuration cannot be achieved with 
counterpropagating waves as in [13 for which the change 
Am = 0, ±2 [22 . Note that it is important that the de- 
tuning is not large compared to the fine structure split- 
ting between the P1/2 and P3/2 manifolds. Otherwise 
the electronic spin S loses its significance, and one is left 
with an effective = ^ Jg = 1 transition where no co- 
herent level superposition is generated by the atom-light 
interaction. 

For the sake of simplicity, we first treat the case of 
a three-level A system. The transposition to a realistic 
internal level structure will be discussed later on. The 
two ground states \gi) and \g2) are coupled to the ex- 
cited state by laser fields with spatially varying Rabi fre- 
quencies 1^1 (r) and ^^2(^), respectively [Fig. [TJa)]. The 
Hamiltonian for an atom in the light field at any point r 
reads 

H=^ + V{r) + Hj,Ur)- (3) 
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FIG. 2: Planar scheme to create a gauge field for atoms with 
two crossed far-detuned laser beams which are displaced by a 
distance iba from the centre position of the atomic cloud. 



The atom is confined in a two-dimensional trapping po- 
tential V{r) = {uolx'^ ^ujly'^)/2M. The part i^AL of 
the Hamiltonian acts on the internal degrees of freedom 
only and describes the atom-laser coupling. Using the 
rotating wave approximation, it can be written as 

/ Vti/2\ 
HAi. = n\ 1^2/2 , (4) 
\^\/2 Vtl/2 -S J 

where S = uoi, — uja '^^ the detuning of the laser frequency 
cjl with respect to the atomic resonance frequency coa. 
We suppose that the atom is prepared in the eigenstate 

IV^oW) =cos(^/2)|^i)+e-^^sin(^/2)|^2) (5) 

of Hal, where we have set cosO = — |172p)/r^^, 

e^^ sin l9 = -2^]^^]2/^^ and Q = + |l]2p)^/^ Here 

is a non-coupled (dark) state, separated from the 
next eigenstate by an energy given hy s = hQ^ /4S in the 
low intensity limit Q <C \S\. The general expressions ([2| 
for the gauge potentials yield in this case 

A{r) = ^(l-cos0)V0, (6) 

U{r) = ^[{V9f + sm'9{Vcj^f]. (7) 

We now introduce a simple implementation scheme to 
generate a gauge field with two far-detuned Gaussian 
laser beams. As the beams will in practice need to be (T_ 
and TT polarised (or equivalently cr+ and tt), we choose a 
crossed-beam configuration. Our beams are displaced by 
a distance ±a (of the order of the beam waists) from the 
X- and the ?/-axis, respectively, as shown in Fig. [2j in or- 
der to obtain a non-vanishing magnetic field at the origin. 
The Rabi frequencies thus read l^i(r) = Qoexjp[ikx — 
{y - a^/w'^] and Q^2(r) = l]oexp[i% - (x + a^/w'^], 
where w denotes the waist of the beams and k the wave 
vector. Expanding the elements of the matrix Q up 
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to second order around the origin and plugging the re- 
sults into Eq. (|6|, we obtain the vector potential A{r) = 
hk{ey — Cx) [1/2 — a{x + y)/w'^] , yielding the effective 
magnetic field to lowest order: 



B{r) = V X A = — —e^ 



(8) 



An important quantity characterising the effect of the 
magnetic field is the cyclotron frequency uOc — q\'B\lM^ 
with = 1 in the present case. For our configuration we 
get = 2fika/Mw'^ close to the origin. For comparison, 
in a gas rotating with angular frequency l^rot one has 
UOc = 2r^rot- In the situation where several vortices are 
present in the ground state, the surface density of vortices 
is related to uOc by [3] 
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The typical radius of the atom cloud can be chosen such 
that R ^ a ^ w so that the number of vortices inside 
the cloud is nR^py ~ ka. Owing to the fast phase vari- 
ation ~ /c ^ 1/a of the light field we may thus expect a 
significant number of vortices to be present in a steady 
state cloud. This number is ultimately limited by the 
higher-order terms in the expansion of the gauge field. If 
the displacement a is chosen too large compared to the 
waist the magnetic field will be inhomogeneous and 
the scalar potential anharmonic. 

To obtain the scalar potential, we note that the first 
term in expression ^ is small compared to the second 
one ~ (V(/))^ = 2k^ and can therefore be neglected. An 
expansion around the origin up to second order then leads 
to 
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which reduces the trapping frequency along the axis x = 
y. This anisotropy can be compensated for by adapting 
the initial trapping potential V. 

Up to now our scheme (like other previous proposals) 
relies on the assumption that the atomic motion is slow 
enough so that the atoms follow adiabatically the state 
IV^o(^))- Although reasonable, this approximation could 
be questioned in the vicinity of a vortex core, where the 
velocity field can become arbitrarily large. To settle this 
point we have numerically determined the ground state 
of a trapped gas of interacting bosonic atoms with a A- 
type internal structure irradiated by lasers in the con- 
figuration of Fig. |2j The gas is assumed to move only 
in the xy plane, and the third motional degree of free- 
dom is assumed to be frozen to the ground state of a 
strongly confining potential along the z direction (fre- 
quency cOz)- Spontaneous emission processes are included 
by adding the complex term —ihT to the energy of the 
excited state. For simplicity we take the same s-wave 
scattering length for all internal states and hence dis- 
regard spin-spin interactions which could lead to spin ex- 
change collisions. This approximation is reasonable for 



the hyperfine ground states of ^''Rb but not necessarily 
for other atomic species. For the excited state it should in 
general be unproblematic as the atoms spend their time 
primarily in the ground state manifold. The interactions 
are treated at the mean-field level, and the ground state 
of the system is obtained using an evolution in imagi- 
nary time of the three coupled GPEs associated with the 
three internal atomic states. A typical result is shown in 
Fig. [3]^ a). We find that the atoms accumulate as expected 
in the non-coupled internal state |'0o(^)) while the pop- 
ulation of states other than is evanescent (fraction 
below 2 X 10~^ for the parameters of Fig. |3|. Further- 
more, the equilibrium COM distribution indeed exhibits 
a vortex lattice whose vortex density at the centre is in 
good agreement with the prediction of Eq. (|9| . The slight 
asymmetry of the cloud shape stems from non-harmonic 
contributions to the scalar potential 

Vortex nucleation following the sudden application of 
a gauge field usually occurs through a turbulent phase 
which lasts a fraction of a second, followed by a relaxation 
to the ground state via a dissipative process [191 EQI |2T] . 
Another possible route could be a smooth evolution from 
a no-vortex state to a state with vortices as the gauge 
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FIG. 3: Density distributions of a two-dimensional (2D) 
gas of bosonic atoms with a A level scheme and rubidium- 
like parameters: resonance wavelength 795 nm, M — 1.45 x 
10~^^ kg, natural linewidth of the excited state r/27r = 
6 MHz. The harmonic trapping potential in the xy plane, 
including the scalar potential [/, is isotropic with a trap- 
ping frequency cj/27r = 40 Hz and a ground state extension 
ttho = (Ti/MijoY^'^ — 1.7/im. Atomic interactions are chosen 
such that G = NasiSnMiv^/ny^'^ = 800. The laser configu- 
ration is sketched in Fig. [2] with a — w = 25 aho, ^o/F = 280, 
and ^/r = 4 X 10^. (a) Steady state distribution, (b)-(d) Dis- 
tributions obtained after ramp durations of the second laser 
beam of (b) 50cj"\ (c) 200 c^"^ and (d) 1000 cj"^ 
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field is slowly turned on. Usually this process is for- 
bidden, at least at the mean field level, because of the 
different parities of the initial and final states. However, 
the gauge field generated in the present scheme does not 
have any definite parity property and such a smooth evo- 
lution should in principle be possible. To test whether 
it can be implemented in practice, we assume that one 
laser beam is switched on first in order to lift the degen- 
eracy between the two ground state levels. Once it has 
reached its full intensity, we ramp up linearly the inten- 
sity of the second laser beam in an adjustable time T. 
Starting with a trapped gas in its ground state in the 
presence of one beam (a Bose-Einstein condensate with- 
out vortices), we propagate the wave function in real time 
using the GPEs. The resulting density distributions are 
shown in Fig.jsj^bcd). For short ramp times T ^ 200 cj"^ 
many vortices are situated around the border of the 
cloud but they have barely moved towards the centre [see 
Figs.jsj^b) andjsjc)]. For slower ramps [Fig. [sjd)] the fi- 
nal state looks closer to the expected ground state shown 
in Fig.jsja), although the convergence is not perfect. We 
conclude that for these parameters the time scale for a 
smooth turn-on of the geometric potential must be longer 
than 1000 4 s for cc;/27r = 40 Hz), which seems too 
long for a practical implementation. Probably vortex nu- 
cleation in these geometric potentials will be more easily 
achieved by following the same route as for a rotating 
trap, i.e., by suddenly applying the gauge field to a con- 
densate at rest and by letting the cloud evolve to its new 
steady state through a transient turbulent phase [19-21]. 

The analytical derivation presented above for a three- 
level system can be straightforwardly extended to the 
case of a realistic atomic transition. Although the re- 
sults remain qualitatively unchanged, it is instructive to 
identify some minor differences. Firstly, unlike for a A 
system, the eigenstates of the atom-laser coupling are 
not "dark" in the sense that they all contain a non-zero 
admixture of the excited level. However, the residual 



photon scattering rate can be made very small. The pa- 
rameters given in the caption of Fig.|3]lead to a scattering 
rate of 0.15 photons/s per atom, which is negligible on 
the time scale of a typical experiment. Secondly, we note 
that the level scheme of Fig. [ijb) leads to some modifi- 
cations of the geometric gauge field with respect to the A 
system. After optimisation of the relative intensities of 
the cr_- and 7r-polarised laser beams, one obtains that the 
effective magnetic field at the origin is multiplied by the 
factor 8J^/(3a/3) with respect to the result For the 
= 1 ground state of ^^Rb, for example, this leads to a 
50 % increase of the vortex density at the trap centre. 

In conclusion, we have proposed in this paper a scheme 
to create a gauge field in a cold atomic gas with negligible 
heating due to photon scattering. We have also validated 
the adiabatic approximation from which this gauge field 
emerges by solving the set of coupled Gross-Pitaevskii 
equations associated with the atomic internal levels. In 
particular, the ground state of a trapped 2D Bose gas 
indeed contains several regularly arranged vortices, with 
the density predicted from the gauge field in the adia- 
batic approximation. The feasibility of our simple scheme 
should render possible an experimental observation of a 
light-induced magnetic field. 
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